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Abstract 

Recently it has been shown that calculations based on a hadronic rescattering model 
agree rather well with experimental results from the first RHIC run. Because of the 
large particle densities intrinsically present at the early time steps of Monte Carlo 
calculations attempting to model RHIC collisions undesirable artifacts resulting in 
non-causality may be present. These effects may compromise the results of such cal- 
culations. The subdivision method, which can remove such artifacts, has been used 
to test the present rescattering model calculations. It is shown that no appreciable 
changes are seen in the present calculations in using the subdivision method, thus 
strengthing the confidence in the results of this rescattering model for RHIC. 
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1 Introduction 



Recent hadronic rescattering model calculations have been shown to agree 
reasonably well with experimental results obtained in the first physics run for 
the Relativistic Heavy Ion Collider (RHIC) with y/s = 130 GeV Au + Au 
collisions [1] [2] . More specifically, these calculations provide good representa- 
tions of the data for 1) the particle mass dependence of the px distribution 
slope parameters (i.e. radial flow), 2) the pr, particle mass, and pseudorapid- 
ity dependences of the elliptic flow, and 3) the px and centrality dependences 
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of the pion Hanbury-Brown-Twiss (HBT) measurements. This is a remarkable 
accomplishment for a single model. 

Many hadrons are initially produced in a relatively small volume in RHIC- 
type collisions, particularly at early times in the interaction. It is a challenging 
task for a Monte Carlo calculation of this type to deal accurately with binary 
collisions between particles in such a large particle density environment with- 
out introducing numerical artifacts which may affect the results. One such 
undesirable artifact that can occur when the interaction range between two 
particles is much greater than the scattering mean-free-path is non-causality 
of the collisions [3] [4]. The method of subdivision can be used to minimize 
these artifacts associated with high particle density [5]. 

The goal of this work is to apply the method of subdivision to the present 
RHIC-energy rescattering calculations to test whether the results from these 
calculations without subdivision are valid. Calculations will be compared with 
and without subdivision for the pion observables 1) (l/pr)dN/dpr, 2) elliptic 
flow vs. pt and r], and 3) HBT vs. pt, at a collision impact parameter of 8 fm. 
Section II will describe the calculational methods used and Section III with 
give the results of the study. 



2 Calculational Methods 



2.1 Hadronic rescattering calculation 

A brief description of the rescattering model calculational method is given 
below. The method used is similar to that used in previous calculations for 
lower CERN Super Proton Synchrotron (SPS) energies [6]. Rescattering is 
simulated with a semi-classical Monte Carlo calculation which assumes strong 
binary collisions between hadrons. The Monte Carlo calculation is carried out 
in three stages: 1) initialization and hadronization, 2) rescattering and freeze 
out, and 3) calculation of experimental observables. Relativistic kinematics is 
used throughout. All calculations are made to simulate RHIC-energy Au+Au 
collisions in order to compare with the results of the y/s = 130 GeV RHIC 
data. 

The hadronization model employs simple parameterizations to describe the 
initial momenta and space-time of the hadrons similar to that used by Her- 
rmann and Bertsch [7]. The initial momenta are assumed to follow a thermal 
transverse (perpendicular to the beam direction) momentum distribution for 
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all particles, 



(l/m T )dN/dm T = Cm T /[exp (m T /T) ± 1] 



(1) 



where itit = Vpt 2 + ?^o 2 is the transverse mass, px is the transverse momen- 
tum, m is the particle rest mass, C is a normalization constant, and T is 
the initial "temperature parameter" of the system, and a gaussian rapidity 
distribution for mesons, 



where y = 0.5 In [(E + p z ) j[E — p z )] is the rapidity, E is the particle energy, 
p z is the longitudinal (along the beam direction) momentum, D is a normal- 
ization constant, yo is the central rapidity value (mid-rapidity), and a y is the 
rapidity width. Two rapidity distributions for baryons have been tried: 1) flat 
and then falling off near beam rapidity and 2) peaked at central rapidity and 
falling off until beam rapidity. Both baryon distributions give about the same 
results. The initial space-time of the hadrons for b = fm (i.e. zero impact 
parameter or central collisions) is parameterized as having cylindrical symme- 
try with respect to the beam axis. The transverse particle density dependence 
is assumed to be that of a projected uniform sphere of radius equal to the 
projectile radius, R (R = t-qA 1 ^ 3 , where r = 1.12 fm and A is the atomic 
mass number of the projectile). For b > (non-central collisions) the trans- 
verse particle density is that of overlapping projected spheres whose centers 
are separated by a distance b. The particle multiplicities for b > are scaled 
from the 6 = values by the ratio of the overlap volume to the volume of the 
projectile. The longitudinal particle hadronization position (zh a d) and time 
(thad) are determined by the relativistic equations [8], 

Zhad = r had sinh y; t had = r had cosh y (3) 

where y is the particle rapidity and r had is the hadronization proper time. 
Thus, apart from particle multiplicities, the hadronization model has three 
free parameters to extract from experiment: <J y , T and Th ad . The hadrons 
included in the calculation are pions, kaons, nucleons and lambdas (n, K, N, 
and A), and the p, u, 77, 77', 0, A, and K* resonances. For simplicity, the 
calculation is isospin averaged (e.g. no distinction is made among a n + , it , 
and 7r~). Resonances are present at hadronization and also can be produced as 
a result of rescattering. Initial resonance multiplicity fractions are taken from 
Herrmann and Bertsch [7], who extracted results from the HELIOS experiment 
[9]. The initial resonance fractions used in the present calculations are: rj/n = 
0.05, p/ir = 0.1, p/oj = 3, 0/(p + u) = 0.12, rj'/rj = K* /u = 1 and, for 
simplicity, A/N = 0. 




(2) 
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The second stage in the calculation is rescattering which finishes with the 
freeze out and decay of all particles. Starting from the initial stage (t = 
fm/c), the positions of all particles are allowed to evolve in time in small 
time steps (dt = 0.1 fm/c) according to their initial momenta. At each time 
step each particle is checked to see a) if it decays, and b) if it is sufficiently 
close to another particle to scatter with it. Isosp in- averaged s-wave and p- 
wave cross sections for meson scattering are obtained from Prakash et al. [10]. 
The calculation is carried out to 100 fm/c, although most of the rescattering 
finishes by about 30 fm/c. The rescattering calculation is described in more 
detail elsewhere [6]. 

Calculations are carried out assuming initial parameter values and particle 
multiplicities for each type of particle. In the last stage of the calculation, the 
freeze-out and decay momenta and space-times are used to produce observ- 
ables such as pion, kaon, and nucleon multiplicities and transverse momentum 
and rapidity distributions. The values of the initial parameters of the calcu- 
lation and multiplicities are constrained to give observables which agree with 
available measured hadronic observables. As a cross-check on this, the total ki- 
netic energy from the calculation is determined and compared with the RHIC 
center of mass energy of y/s = 130 GeV to see that they are in reasonable 
agreement. Particle multiplicities were estimated from the charged hadron 
multiplicity measurements of the RHIC PHOBOS experiment [11]. Calcula- 
tions were carried out using isospin-summed events containing at freezeout for 
central collisions (6 = fm) about 5000 pions, 500 kaons, and 650 nucleons 
(A's were decayed). The hadronization model parameters used were T = 300 
MeV, a y =2A, and 17^=1 fm/c. It is interesing to note that the same value of 
T had w as required in a previous rescattering calculation to successfully describe 
results from SPS Pb+Pb collisions [6]. 



2.2 Subdivision method 

The method of subdivision is based on the invariance of Monte Carlo particle- 
scattering calculations for a simultaneous decrease of the scattering cross sec- 
tions by some factor, /, and increase of the particle density by Z[3]. As I becomes 
sufficiently large, non-causal artifacts become insignificant. The present rescat- 
tering calculation will be tested comparing pion observables for no-subdivision, 
i.e. 1 = 1, with subdivision of I = 5. A subdivision of / = 5 was chosen since 
when applied to other Monte Carlo scattering calculations this value has been 
shown to produce large changes in observables when non-causal artifacts are 
present [3] [4]. Since the particle density increase is accomplished by increasing 
the particle number by a factor /, the computer CPU time taken per event 
increases by a factor I 2 , effectively "slowing down" a subdivision study by 
a factor I compared with not using subdivision. For the present study, 1640 
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events and 124 events were generated for the I = 1 and I = 5 samples, re- 
spectively (the statistical value of the I = 5 events being 5 times greater per 
event). The CPU time taken to generate the I = 5 sample of events was 420 
CPU-hours on a 1 GHz PC processor. 



3 Results 

Figures 1-4 show comparisons of I = 1 with I = 5 hadron rescattering model 
calculations for pion observables. Descriptions of how the observables are ex- 
tracted from the rescattering calculation are given elsewhere[2]. All calcula- 
tions are carried out for an impact parameter of 8 fm to simulate a medium 
non-central collision which should result in significant elliptic flow for the pur- 
poses of the present test. Statistical errors are shown either as error bars or 
are of the order of the marker size when error bars are not shown. 

Figure 1 compares the pt distributions for a rapidity cut of —2 < y < 2. As 
seen, the / = 1 and / = 5 calculations agree within statistical errors, although 
there is a slight decreasing trend for / = 5 compared with I = 1 for p T > 3 
GeV/c. 

Figures 2 and 3 compare elliptic flow (V2). Figure 2 shows V 2 vs. pt for a 
rapidity cut of — 2 < y < 2 and Figure 3 shows V 2 vs. i] for a px cut of 
< pt < 3 GeV/c. For both cases, the / = 1 and I = 5 calculations are seen 
to agree within statistical errors, reproducing the "flattening" in V 2 at large 
p T and the "peaked" dependence of V 2 on rj. 

Finally, Figure 4 compares HBT results vs. pt for a rapidity cut of — 2 < y < 2. 
Each result is calculated in a pt bin of 0.2 GeV/c. A three-dimensional fit is 
made to extract the pion source parameters RTside, Rrout, Riong, and A. As 
seen in this figure, the / = 1 and I = 5 calculations agree within statistical 
errors for all parameters at all pt bins, reproducing the decreasing trend of 
the i?-parameters and increasing trend of the A-parameter with increasing pt- 

It can be concluded from Figures 1-4 that non-causal artifacts of the type 
which can be removed by the subdivision method are not present in the 
hadronic rescattering code used in this study. It is possible to speculate on 
why this is the case. There are three main features of the code which might 
tend to reduce these artifacts: 1) individual particles are allowed to scatter 
only once per time step, 2) a "scattering time" of two time steps is defined 
during which particles that have scattered are not allowed to rescatter, and 
3) once two particles have scattered with each other, they are not allowed to 
scatter with each other again in the calculation. The present study was per- 
formed for non-central collisions (6 = 8 fm), but it is expected that the same 
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Fig. 1. Pion distributions for I = 1 and 1 = 5. 

conclusion, i.e. that subdivision does not significantly effect the results of the 
calculation, would also be obtained for central collisions. This is because the 
particle density is not very different for central collisions in the model than 
for mid-peripheral collisions since the particle multiplicities are scaled by the 
overlap volume for b > 0, as mentioned earlier. 

In summary, the subdivision method has been used to test the validity of 
the present rescattering model calculations. It is found that no appreciable 
changes are seen comparing I = 1 and I = 5 calculations, thus strengthing the 
confidence in the results presented previously from this rescattering model for 
RHIC. 
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Fig. 2. Pion elliptic flow vs. px for Z = 1 and / = 5. 
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Fig. 3. Pion Elliptic flow vs. pseudorapidity (77) for I = 1 and Z = 5. 
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Fig. 4. Pion HBT vs. px for I = 1 and / = 5. The ordinate scale for A is shown to 
the right. 
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